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Jets are expected to play a prominent role in the ongoing efforts to characterize the 
hot and dense QCD medium created in ultrarelativistic heavy ion collisions. The success 
I of this program depends crucially on the existence of a full theoretical account of the 

• • dynamical effects of the medium on the jets that develop within it. By focussing on 

^ the discussion of the essential ingredients underlying such a theoretical formulation, we 

aim to set the appropriate context in which current and future developments can be 
rS understood. 

^ Keywords: Jet physics; jet quenching; heavy-ion collisions. 

PACS numbers: 12.38.-t, 24.85. +p, 25.75.-q 



1. Introduction 

Experimental measurements at RHICJSEI and the LHCQ have established extensive 
and unequivocal evidence supporting the expectation that a deconfined state of 
matter, commonly referred to as a quark-gluon plasma (QGP), is created in the 
aftermath of ultrarelativistic heavy-ion collisions. This hot and dense state provides 
a unique setting in which dynamic and collective properties of QCD can be studied 
under conditions resembling those prevalent in the primordial Universe. 

'Corresponding author. 

^For a recent review see, e.g., Ref. pjand references therein. 
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Rather than relying on a single 'smoking gun' signature, the study of the proper- 
ties of the QGP resorts to a wealth of complementary measurements which include: 

• the behavior of low-pjj^hadrons, the bulk outcome of heavy- ion collisions, 
from which collective properties in both the initial conditions (the satura- 
tion of the partonic densities in the colliding nuclei) and the final state (the 
development of hydrodynamic flows) can be inferred; 

• the modification of the formation probability of quarkonia states due to the 
presence of a QGP; 

• the modifications, the focus of this review, effected by the QGP on the 
high-pT debris produced in the collision. 



As realized long ago by Bjorken,^ ^ 'hot spots' — or rather hard partonic colli- 
sions — occur unaffected by the plasma due to the short timescales involved. Thus, 
the high-pT products of these processes can be reliably calculated by perturbative 
methods (see Sec. [2|. The effect of the plasma on these 'hard probes', commonly 
referred to as 'jet quenching', and their potential use to extract medium properties 
have been actively studied for the best part of the last three decades (see Refs.jS - 11 
for early pioneering works and Refs. 12- 15 for recent reviews). 

The advent of the LHC heavy ion program — with a large increase of collision 
energy (2.76 TeV/c per nucleon pair as compared to the 200 GeV/c available at 
RHIC) and the excellent detector capabilities of ALICE, ATLAS and CMS — have 
given access not only to a much extended kinematic range for those observables 
previously explored at RHIC, but also to a whole range of novel measurements in- 
volving jets reconstructed reliably and systematically from the midst of the large 
and fluctuating backgrounds created in nuclear collisions. The availability of these 
measurements undoubtedly calls for significant refinement of the underlying theo- 
retical descriptions. 

Within the limited space available for this survey we, rather than providing an 
exhaustive overview of developments and results obtained over the last decades (for 
this purpose we refer the interested reader to the references above), discuss, arguably 
in a new light, the essential ingredients necessary for quantifying the modifications 
experienced by jets in ultrarelativistic heavy ion collisions. Early efforts in this 



direction can be found in Refs. 16 17 



The aim of the 'hard probes' program in heavy-ion collisions — to extract infor- 
mation about the properties of the soft and deconfined QGP — is pursued via the 
identification of deviations from a well-calibrated baseline established in the absence 
of a medium, i.e. in proton-proton collisions. This generic strategy can be illustrated 
with its application to the case of the nuclear modification factor Raa which eval- 
uates the deviation of single-particle inclusive spectra away from the baseline and 



''Throughout, will denote transverse momentum in the detector frame. 
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is given by 
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where X represents the species of the probe and Nab and Npp are, respectively, 
the inclusive yields in nucleus- nucleus [A + B) and proton-proton {p + p) collisions 
{dVL is simply the relevant phase space measure). The ratio is scaled by the (av- 
erage) nuclear overlap function (Tab), which corresponds roughly to the number 
of elementary nucleon-nucleon collisions expected for any given centrality class of 
the nucleus-nucleus collision and is usually extracted from a probabilistic Glauber 
modelEllEHl The ratio is defined such that Rab = 1 when no modification is present. 
The nuclear modification factor is, in general, a function of the collision centrality 
(equivalently the impact factor 6), transverse momentum, etc. of the measured par- 
ticle X. 

The same procedure can be used to devise more general nuclear modification 
factors that probe two- (and higher-) particle correlations or fully reconstructed 
jets. Clearly the range of relevant observables is not exhausted by such ratios. In 
particular, a significant part of current emerging picture of how jets are modified 
by the medium has been obtained from observables — the energy imbalance in 
dijet pairs ,12211211 between isolated photons and their recoiling jet P^IMI and between a 
high-pT trigger hadron and its recoiling jetp^the geometricaP^l and momenturrPmUl 
fragmentation patterns of jets — for which information is more readily conveyed 
through the separate consideration of the measurements in heavy ion collisions and 
in the baseline than by their ratio. 

To attribute the observed modifications, in jet and jet-like hadronic observables, 
to the presence of a QGP it is essential to account for analogous effects arising 
from other sources: the nuclear modification of the partonic distributions in the 
colliding nuclei (an initial state effect); and modifications already present in low- 
energy nucleus-nucleus and in proton-nucleus collisions, both cases where no QGP 
is formed. For the high-py observables on which this review is focused, the nuclear 
parton distribution functions are well constrained over the relevant kinematical 
range. The effects observed at low energy and in proton- nucleus collisions are due 
to purely hadronic, essentially non-perturbative, scattering mechanisms. As such, 
many of their characteristics — species (meson/baryon), transverse momentum and 
rapidity dependences — are not understood from first principles but are well con- 
strained by existing experimental dataESEIl Generically, these 'cold nuclear matter' 
effects are modest and expected to vanish at high-p^- Indeed, in pPb collisions at 
LHC energies (i/s = 5.02 TeV per nucleon pair) the nuclear modification factor is 
close to unity for pr > ^ GeV/cW^ 

Throughout this review we address the interaction of the hard probes with the 
QGP at the partonic level. Such a perturbative treatment is motivated by the large- 
ness of the typical probe momentum scales as compared to the ones characterising 
the medium. This separation of momentum scales implies a separation between the 
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timescales related to the probe and the long-distance features of the medium. In 
other words, the probe (or specific properties of the probe) is only sensitive to those 
medium characteristics involving similar timescales. 

One should note, however, that conceptually orthogonal approaches where the 
probe-medium coupling is taken as large, and necessarily non-perturbative, have 
been actively pursued in the context of the AdS/CFT correspondence (for a recent 
and comprehensive review, see Ref. 33 1 . 

In summary, this review addresses what we believe are the essential theoretical 
ingredients for the description of the development of jets in the presence of hot, dense 
and colorecQ medium. We restrict ourselves to the discussion of jets initiated by 
massless partons, that interact perturbatively with the medium, and whose medium 
modification is dominated by radiative processes. 

The many approaches devised to describe medium induced radiation (BDMPS- 
ZjSHm GLV^ HTplEll AMYEsEZI and their extensions) use very similar assump- 
tions about the main features of the process under consideration but differ in the 
detailed implementation of medium effects and in approximations allowing for an- 
alytic treatment. Although these differences play a significant role when comparing 
to experimental data,!^ we refrain from presenting them as competing implemen- 
tations, but rather attempt to find a common ground from which progress can be 
made and areas where increased rigor is necessary. 

Vacuum jet physics is discussed in Sec. [2] with particular emphasis on the aspects 
which will be object of the medium induced modifications addressed in this review. 
In Sec. [3] we describe the medium we take as input for the jet modifications detailed 
in later sections: medium induced radiation in Sec. [4] and the medium effect on the 
color properties of a jet in Sec. [5j We conclude in Sec. [6] with a brief assessment of 
the phenomenological status of jet quenching. 



2. Jets in vacuum: baseline and template 

Starting with the factorization property of QCD cross sections, we discuss 
in brief the essential features of jet physics in vacuum: the collinear nature 
of gluon radiation and color coherence. 

The essential baseline, with respect to which the modifications of jet observables 
due to the creation of a QGP can be identified, is provided by the description of 
the same observables in proton-proton collisions. In this case, where jets develop in 
vacuum, the current level of understanding is of enviable maturity and precision. 
Jet properties are reliably computed both analytically and in Monte Carlo event 
generators. 

To a great extent, the factorizability of short- and long-distance processes in 
QCD underlies these successes. Such factorizability is evident when the production 

•^We emphasize that since a hot QGP is a poor color conductor,!^ the partons that propagate 
through it will exchange color charge at a rapid rate. 
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of a final state X (parton, hadron, jet) is written in the form 

®Dk^xiz,Q^). (2) 

Here, the PDF //^i (xi, Q^) ^fh2(^^^^Q2^^ accounts for the probability to find a 
parton of species i (j) within the incoming proton hi (/12) carrying a fraction Xi 
{X2) of the total longitudinal momentum and of virtuality (the hard scale set by 
the partonic process). PDFs are universal (process independent) non-pertubative 
objects with scale evolution driven pertubatively by the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equationsP^HUlxhe relevant initial conditions are 
determined via global fits to data and, for the kinematics relevant for jet production, 
are reasonably well constrained. 

The hard scattering of partons i and j, described by the perturbative cross sec- 
tion cr*^"*''^, results in a partonic system (a pair of back-to-back partons at leading 
perturbative order) containing the parton k we focus on and which will relax its 
initial virtuality through QCD branching down to a scale Qq below which per- 
turbation theory ceases to be applicable and hadronization of the fragments takes 
place. 

The probability for the (time-like) parton k (quark or gluon) to branch into 
partons / and m carrying respectively fractions z and 1 — z of its momentum is 
given by 

dw''^'+^" = ^^dzP,,{z), (3) 
47r fe^ 

where the Pik{z) are the Altarelli-ParisP21 splitting functions which also describe 
the decay of a negative virtuality, i.e., space-like evolution of PDFs, and fcj^ is the 
transverse momentum generated in the branching process. 

The well known fact that QCD branching is dominated by the emission of soft 
gluons at small angles is explicitly seen in Eq. (|3|: the logarithmic enhancement of 
small angle radiation from / (Pk± k'J^ cx logQ'^/Q^; and the soft enhancement from 
the integration over z of the contributions in the splitting functions where a gluon 
carrying a small momentum fraction is produced (cx J dz/{l — z) with z ^ 1 and 
cx J dz/ z with z — > 0). 

The space-time structure of multiple branchings in the logarithmic regions is 
characterized by a strong ordering of the typical time scales of successive branchings, 
t{ ~ z(l — z)E/k\. Unlike the space-like evolution, where the strong ordering in 
transverse momenta fcj^ accounts fully for the resummation of large logarithmic 
enhancements, a^. log Q^/Q^ ^ 1, the time-like evolution is determined by a strict 
angular ordering of successive branchings due to color coherence effects: 

Take the outcome of a branching, a pair of partons with angular separation 9. 
The radiation of a further parton of energy w at angle 9 from its emitter, in other 
words a successive branching, takes a finite time which can be estimated from the 
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uncertainty principle as if ~ {uj9'^)~^. Since w0 = /c_l = AJ^, with A_l the wavelength 
of the radiated parton, the formation time can be rewritten in the form if X±/9. 
Within this time, the original partons separate a distance r± ^ 9 1{ ^ X±9/9. Hence, 
for large angle emissions 9 ^ 9, r±^ < Aj^, the radiated quantum cannot resolve the 
structure of the initial parton pair and probes only its total color charge. Conversely, 
small angle radiation occurs from the larger color charge of the resolved partons in 
the pair and is consequently favoredEHIHI 

The resumation of logs leads to a set of coupled evolution equations of the 
Modified Leading-Logarithmic Approximation (MLLA) fornP^tSU 

^A(.,<3)^i:/'!|2#^P„(=)D,W«.«<3), (4) 

where k± = z(l — z)Q, for the fragmentation functions Dj, where now Pji{z) stands 
for the regularized splitting functions in the limit z — 1 after the inclusion of virtual 
corrections. The leading effects of color coherence are accounted for by the shifted 
scale zQ in Dj on the right-hand side of Eq. Q which effectively implements the 
angular ordered pattern argued above by restricting the angular range available 
for each splitting. We recall that in standard DGLAP evolution ,1^21111] where color 
coherence is not accounted for, the scale in Dj is simply Q. While unimportant 
for most values of x, this modification is significant at small x and essential for 
the description of the experimentally observed humpbacked plateau in the inclusive 
energy spectrum of particles in a jet. 

Together with the dominance of soft radiation underlying Eq. Q, the angular 
pattern dictated by color coherence defines the space-time structure of jet develop- 
ment: branchings occur successively with decreasing emission angles, the so-called 
angular ordering of radiation, and are increasingly softer. 

So far we have discussed the fragmentation of the initial hard parton k into a 
partonic system. It is not these 'final state' partons that are experimentally ob- 
served, but rather the hadrons they give rise to. The simple hypothesis of Local 
Parton-Hadron Duality (LPHD^P^ which lies in the assumption that the conversion 
of partons to hadrons affects inclusive observables only via an overall multiplica- 
tive factor and the non-pertubative scale Qo at which perturbative evolution is 
interrupted, has been very successful; see alscEZEH for calculations of intra-jet fcj^- 
distributions invoking LPHD. However, for more refined analysis one has to rely on 
hadronization models to deal with different particles species or for unfolding parton 
distributions from hadron spectra. 

While the, intrinsically non-perturbative, physics of hadronization is not under- 
stood from first principles, it can be effectively encoded by generalizing the fragmen- 
tation fmictions D^. Fragmentation functions Di^x{x,Q) describe the probability 
distribution for a generic final state X (hadron, jet, ...) to result from the branch- 
ings and subsequent hadronization of fragments of a parton i of which it carries a 
fraction x of longitudinal momentum. Their evolution is driven by Eq. m and their 
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non-perturbative input is constrained by global fits to data in a manner analogous 
to the PDFs. 

Alternatively, as is the case Monte-Carlo event generators, the partonic branch- 
ing can be performed down to a scale Qo ^ 1 GeV at which an effective hadroniza- 
tion dynamical procedure is invoked. In these implementations, the partonic frag- 
ments are grouped into color neutral structures (Lund strings,!^ clusterJ^ which 
dynamically decay into the final state hadrons. 

When addressing observables involving reconstructed jets, a strict definition — 
an algorithm specifying how to group the fragments completed by a set of param- 
eters (e.g., the jet radius) — of what the jet is must be given. Since, in general, 
different jet definitions result in different jets, comparison between theory calcula- 
tions and data are only meaningful when the same definition is used. The rationale 
behind the various existing jet definitions, their applicability and robustness are 
discussed at length in Ref . [61] 

3. Probing the medium 

The jets probe the underlying medium, which we proceed to discuss. We com- 
pare different models for the medium and consider thermal effects that arise 
in the plasma. Subsequently, the concrete model realization of the medium 
properties is treated as an input to the calculation of medium effects on jets 
in the following sections. 

In the rest frame of a highly energetic particle the incoming medium is strongly 
Lorentz contracted and nearly translationally invariant. One therefore typically as- 
sumes that the probe will not be sensitive to the longitudinal structure of the plasma 
but only to its static properties. In other words, that the interaction between probe 
and medium is instantaneous. Leaving a discussion about this point to the end of 
this section, we will presently implement this approximation which translates into 
the fact that the momentum exchange is purely transverse. The medium gauge field, 
^med(9) = ^°^mcd('?)0^^^^^ 9 momentum transfer from the medium and 

is a SU(3) matrix in the fundamental representation, can therefore conveniently be 
written in terms of the effective field 



which, in this mixed representation, depends only on the light-cone time and 
the transverse momentum q = {qi,q2). To simplify the expression in the following 
we will identify the symbol t with the light-cone time, t = x'^ . 

In the setup described above, assuming a Gaussian distribution for the medium 
field, the input from an underlying theory of the plasma enters in the simplest case 

"^We will work in the light-cone gauge A'^ = where A~ is the only relevant component of the 
medium field. 

"^Light-cone coordinates are defined as = (x" it x'^)/%/2. 
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as a two-point function correlator, which can be written as 

«;cd(g. O^'cd" (q', t')) = 6^' nit) 5{t - t') (27r)2<5(2) {q - q') 7(9^) , (6) 

where 7(9^) contains the microscopic details of the interaction with the medium 
constituents and n{t) the density of color charges (which could be a function of 
the interaction time for expanding media). Note that the correlator is instanta- 
neous reflecting the assumption about a translationally invariant medium. Neglect- 
ing higher-order correlators in all observables corresponds to treating the medium 
as a set of independent scattering centers. In fact, this correlator, at lowest order 
in the medium coupling scales with the medium length and is leading compared 
to higher-order ones in the limit of large media and g <^1^^ 

This can also be understood in terms of screening phenomena in the plasma. 
Consider the squared Yukawa potential (in momentum space) 

7Gw(g^) = , ^^\. , (7) 



[q^ + to| 



which is screened by the characteristic (Debye) mass to_d. In coordinate space this 
corresponds to the fact that the potential extends up to a characteristic distance 
Tscr ~ '^'d^ ■ This physical setup models the medium as a set of static, randomly 
distributed scattering centers with a mean free path given by A^fp ~ (n(i) CTci)"^, 
where the elastic cross section is simply (Tei oc / (Pq^fiq^)/ (27r)^. This is the so-called 
Gyulassy-Wang (GW) model of the QGPlHlllllIll Assuming that rscr < A^fp allows 
us to treat the scattering centers as independent and justifies the simplifications 
above. In the opposite case the probe can, in principle, be sensitive to higher-order 
correlators which capture collective behaviors of the plasmaj^ 

Gluon fields at finite temperatures generate screening effects as well. These ef- 
fects can be studied by high-temperature effective theories, such as the hard thermal 
loop (HTL) approximatiorHUmi (see, e.g., Refs. 68p9 for reviews). It was found that 



the longitudinal (electric) gluon fields are screened by a dynamically generated De- 
bye screening mass, which relates to the temperature of the medium as m^i ~ gT. 
The (static) magnetic components are, on the other hand, not screened (see, e.g.. 



Ref. 70 for further details). Since the mean free path scales with an additional factor 
of the inverse density, and thus scales as {g'^T)~^, in the weakly coupled regime, 
(7 1, the assumption of independent scattering centers is justified.!^ 

First and foremost, the interaction with the thermal fields induces thermal 



masses of the probe due to the modification of their (static) self-energies, see Eq. ( 23 ) 
below. Besides, the probe can alter its kinematics during propagation by exchanging 
momentum with the thermal medium. Historically, one first calculated the elastic 
rescattering cross sectiorEDEH which gives rise to so-called elastic energy loss-^^Due 
to the inherent collinear nature of radiative processes it was quickly realized that, 
although being suppressed by a power of the coupling, they could contribute at the 



^Non-eikonal corrections, allowing for recoil effects, have recently been calculated in Ref. 65 
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Fig. 1. Diagrams contributing to the dipole scattering rate in Eq. (one must also add the 

complex conjugate diagrams). 

same level as elastic scattering due to phase space enhancement and become domi- 
nant for propagating partons of sufficiently high energy. This was first discussed in 
the context of photon radiation at finite temperatureJMlE^EIlEIl and later extended 
to gluonsElEZlThis gives, in turn, rise to the so-called Landau-Pomeranchuk-Migdal 
(LPM) effectCSEllthe formation time of induced radiation can exceed the mean free 
path giving rise to interference effects between subsequent rescatterings, see Sec. |4.1| 
for a comprehensive discussion. Since the radiative processes scale with a larger 



power of the in-medium path length, see Eq. ( 29 ) and discussion below, compared 
to elastic ones, one usually neglects the latter effects for highly energetic probes and 
large media. ^ While elastic rescattering effects should be incorporated consistently 
for low-pr observables, see also Ref. '78' and comment below, we will not currently 
examine them in more detail|3 

Then, for soft momentum transfers from the medium, |q| ^ T, the potential 
(squared) at leading order in the coupling become^^ 

7HTL(q') = 2, 2 ^ ' (8) 

q2(q2 + 



and scales as 7htl ^ Nq~'^ for |q| T, where the constant is e.g. given in Ref. 84 
Comparing to the static potential, Eq. ([t]), one observes a divergent behavior for 
small \q\. Higher-order corrections in g to Eq. ^ are also knownj^^and lead to an 
even bigger enhancement of the soft sector. Thermal effects are included in several 
theoretical calculationiMtt£Li.2^L§^ of radiative processes in medium, recently also in 
the presence of a finite chemical potential.'^ 

From our discussion so far, the probe will be sensitive to medium characteristics 
through interactions which induce dependence on parameters. The second moment 
of the correlator in Eq. ([6]), historically called g, is a measure of the transverse 
momentum (squared) acquired by the probe per unit length in the elastic scattering 
and, as we will see below, is a highly important quantity for the study of jets in 
medium. We will define it stripped of its relevant color factor, as 

q{t) = (^sn{t) j dq'q^-iiq^), (9) 

SQther effects, such as, e.g., transition radiation,l2ll absorptive efTecti§2I£U and Mach cone cre- 
atioilSSIor Cherenkov radiatiorPS due to supersonic motion in the plasma, can also play a role but 
will not be discussed here. 
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where q* is a ultraviolet (UV) cut-ofF and 7(<7^) is given by the medium-model 
under consideration, cf. Eqs. Q or Q.*^ Anticipating the detailed discussion in 
Sec. |4j the relevance of q can be justified by considering how a dipole of transverse 
size r interacts via two-gluon exchange with the medium, see Fig. [l] The resulting 
scattering rate, stripped of its relevant color factor, is given by 

T2ir, t) = 2gMt) j - e'"") 7(9^) • (10) 

To leading logarithmic order one finds back q by expanding the term in the brackets 
to the first non-trivial order and comparing Eq. ([9]) with the previous equation it 
becomes clear that the relevant cut-off q* ~ I'^'l""'^ in Eq. (lIO).*^ 

For both the static and HTL forms of the potential, Eqs. (7| and ([8|, the general 
form of the scattering rate is r2 oc r^(log l/(TO|,r^) + const.). In the limit r^^ ^ 
rn\) one can neglect the variation of the logarithm which corresponds physically to 
multiple soft scatterings in the medium. In this limit we can employ the so-called 
'harmonic oscillator' approximation which reads 

T^{r,t)c^\mr\ (11) 

and has been historically used to ease the analytical treatment of the expressions. 
The additional logarithmic behavior in the opposite limit arises due to hard mo- 
mentum transfers with the plasma constituents (quasi-particles). This situation is 
in fact highly relevant for highly energetic probes since the UV regulator q* scales 
with a power of its energy j^ ^'^n^^a^'-'. This 'hard' contribution has a strong impact on 
phenomenological observable^^ and should therefore be treated carefully. Keeping 



this in mind, throughout we will mostly employ Eq. ( 11 ), commenting on the 'hard' 
dynamics whenever it is relevant. 

As will be demonstrated explicitly in Sec. |4) for energetic (eikonal or close-to- 
eikonal) particles in the medium, the only parameter — apart from the thermal 
masses — that governs the medium-induced dynamics is q defined in Eq. (|9| . It 
determines the transverse momentum broadening of particles in the medium, see 



Eq. (39), and controls the rate of medium-induced splittings, Eq. ([26|. This parame- 
ter can be calculated from an underlying microscopic theory, e.g., high-temperature 
QCD in the HTL approximation in isotropic!^ or anisotropic!^ background, or on 
the lattice,'^ but it is most usually treated as a phenomenological parameter to be 
extracted from datalHUUll 

Above, cf. Eq. ([9]), g is defined locally and, as such, should be understood as a 
transport coefficient of the medium. In the literature one has also studied it globally 
in the context of the transverse momentum broadening of a probe traversing a 
medium of length L, see also the discussion in SeclIilandEq. (Isi) HHMZl Within 



the approximations discussed above (translationally invariant medium, etc.) these 
two definitions are equivalent, but the global approaches opens up the possibility 

''Since |r|~-' ^ fe^ri see below Eq. l |l4| , we expect q* oc _E^/*. 
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for studying the interplay of medium scale^HZHHl and long-distance behavior in more 
detail. In this case one could also account for the recoil of induced radiation in 
course of the propagatiorP^l in the definition. 

Finally, we note that the assumption of a translationally invariant medium 
should break down for probes with lower energy which subsequently become sensi- 
tive to the longitudinal structure of the medium. This gives rise to energy depletion 
of the probe or so-called drag effects. In the most general case, the correlator in 
Eq. ^ is related to the retarded propagator of the medium fields,!^ and there- 
fore becomes responsive to other moments, i.e. not only involving the transverse 
momentum components, besides q. These novel components are called e and 62 
and can be incorporated in analysis of the propagation of a color probe through 
the medium .lifiSI They provide an additional source of energy loss which is closely 
related to the elastic energy-loss mechanism mentioned above. 

It is worth recalling that the above parametric estimates are valid in a weakly- 
coupled regime, where g <^ \. Going beyond this approximation is very challeng- 
ing from a theoretical point of view.^^^ In a complementary effort, strong-coupling 
techniques based on gauge-gravity duality have recently been used to address the 
problem of energy loss in plasmas; for recent reviews see, e.g., Refs. |33l[l02| and ref- 
erences within. Ultimately, one aims at controlling the remaining dynamics of both 
probe and heavy-ion system such that one can reliably extract information about 
medium characteristics on the microscopic level. An interesting perspective in this 



sense is, e.g., to note the different UV behavior of Eq. (10 1 in weakly and strongly 
coupled plasmas,!^ which arises due to the absence of quasi-particles in the latter 
regime. 



4. Induced radiation in the medium 

Interactions with the deconfined medium induce the radiation of gluons off 
the jet. We discuss how this effect probes the medium via the so-called LPM 
effect and show how this radiation becomes a source of energy loss. Contrary 
to the collinear nature of radiation in vacuum, this component tends to 
propagate to large angles and decoheres in color from the remaining jet 
structure. 

The suppression of high-px particles observed in high-energy heavy-ion collisions 
was early interpreted in the context of radiative energy loss.l^W^ These calculations 
predicted that energetic partons tend to lose their energy via stimulated radiation 
of soft gluons as a result of interacting with the color charges of the dense medium 
formed in the collision. The resulting shift in energy, Pt ^ Pt ~ ^PT) leads to a 
depopulation of high-pj- modes. Our goal is to discuss the state of the art of the 
theory of jet-quenching which is common to most of the present-day calculations, 
its achievements and limitations, and provide the reader with the basics tools and 
references. 
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4.1. The LPM effect in QCD 

Consider an energetic parton of energy E, produced in a hard process very early 
after the collision, typically over a time 1/E, traversing a medium of length L. It will 
undergo multiple final state interactions with the color charges of the dense medium. 
Given that the energy of the parton is much higher that the typical momentum 
acquired due to the multiple scatterings in the medium, i.e., E :S> {p™°'^{L)), the 
parton-medium interactions mainly cause color precession at high rate.!^ These 
color excitations stimulate gluon emissions continuously along the in-medium path. 

As a naive expectation for the emission rate let us assume that all the A'gcatt = 
L/Amfp scattering centers along the path act independently. This implies that the 
spectrum will be an incoherent sum of emissions due to every scattering and reads 

^jindcp 

^ — ^ ~ a^A^scatt • (12) 

auj 

This estimate is of course wrong because it ignores coherence effects due to quantum 
interferences between scattering centers. In fact, medium-induced gluon radiation 
can be understood in terms of medium resolution power of the higher Fock-states 
of the high-energy parton. 

As discussed at length in Sec. [3j the medium is characterized by the parameter 
q which corresponds to the characteristic transverse momentum squared per mean 
free path, so that q ^ nr^/X^fp heuristically. A typical fluctuation of the quark 
state, namely a virtual gluon lives a time At and can explore the plane transverse 
to the projectile over a distance Aa;_L ~ i/k± due to quantum diffusion. During 
this time, if the average transverse correlation length in the medium (qAt)^^^^ 
— which is simply the typical inverse transverse momentum that a particle can 
accumulate during At — is larger than the size of the fluctuation, the fluctuation 
will be transparent to the medium. But as soon as Ax±_ ^ {qAt)~^/^ , which implies 
that fc^ ^ qAt, the medium resolves the system gluon-quark and the gluon is 
freed. Making use of the uncertainty principle. At ^ k~ ^ ui/k\, we then obtain 
an estimate for the characteristic timescale for this process, which we can call the 
gluon radiation time 

tbr-y|, (13) 

which also gives the characteristic medium-induced momentum 

fcg, = (14) 

accumulated during <bi - This gives consistently the size of the gluon at formation 
as Ax±^ ^ ^hr ■ follows that gluons traversing the full length of the medium 
acquire the maximal amount of accumulated transverse momentum, denoted Qs ~ 
y/qL. The scattering centers involved in the radiation process act coherently as 
one. Thus, even if the medium is dense the actual in-medium mean-free path of the 
fluctuation is no longer the elastic one Amfp but rather <br due to color transparency. 
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Therefore the number of effective number of scattering centers is N^s = i/tbr- This 
quahtative discussion allows us to correct the above naive estimate of the rate of 



gluon emissions, by replacing iVscatt N^s in Eq. (12 1, so that 



— w— a^W — , (15) 

which involves the characteristic gluon energy = qL^ ior giuons with tbr = L. 
Thus, coherent medium-induced radiation leads to a suppression of the hard modes 
scaling as as compared to the incoherent case. This form of the spectrum. 



Eq. (151, is expected to hold between two critical frequencies wbh < w < w^, where 
we have defined the so-called Bethe-Heitler frequency ojbh = ^'^mfp- This frequency 
characterizes the giuons that are radiated over times of the order of the elastic mean 
free path — giuons that probe smaller lengthscales are naturally not expected to 
be induced in the medium. Note that this lower cut-off on energy ensures that 
the giuons are emitted in the forward direction.ESSl Indeed, demanding that the 
emission angle is smaller than 7r/2, or fcbr < ^, leads to a; > wbh > <?^^^, where the 
last inequality follows from the condition of independent scatterings. 

This phenomenon is the QCD analog of the Landau-Pomeranchuk- 
MigdadZHEnil effect in QED that leads to the suppression of the soft photon 



frequencies, see e.g. Ref. 105 



4.2. In-medium branching rate 

To see how the physical picture discussed above arises concretely in the calculations 
it is worth to quickly recall the basis for the techniques that have been developed 
to calculate radiative processes in medium. In this section we discuss the leading 
order medium-induced gluon emission. 

We shall assume that energies of the quark and the radiated gluon are much 
larger than their transverse momentum or, equivalently, to the momentum trans- 
ferred to the medium, so that ^ p± and k'^ ^ A;_Lj^This kinematics corresponds 
to the so-called eikonal approximation of the projectile-medium interaction in which 
the field theoretical description becomes equivalent to that of two-dimensional non- 
relativistic quantum mechanics. Also, in the following we shall restrict our discus- 
sion to the purely gluonic case since the generalization to include quark degrees of 
freedom is straightforward. 

Working in the light-cone gauge, A+ = 0, we can show that only transverse gluon 
polarizations propagate, see e.g. Ref. 106 In addition, the eikonal approximation 
implies no spin flip along the trajectory yielding a diagonal propagation of the gluon 
polarization. Thus, in the presence of a two-dimensional background field ^i„cd(^' 
see the left-hand-side of Eq. ([5]), the gluon propagator Q{x,t;y,t') between two 



'The momenta p and k denote the final-state momenta of the projectile and radiated gluon, 
respectively. 
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d 



*.9A„cd(^'*) 



g{x,t-y,t') = i5{t-t')5{x-y). 



(16) 



light-cone times t and t' and transverse positions x and y, respectively, obeys the 
following Schrodinger-like equation 

2p+ 

where the transverse position is conjugate to the transverse momentum p± and the 
light-cone momentum p+ plays the role of a mass and is conserved in the eikonal 
propagation. Note that in the strictly eikonal limit (p"'" ^ p±) the solution for 
Q{x,t;y,t') is a Wilson line along the trajectory of the particle, given by x{t) = 

y + {t — t')p/p'^ , such that 



Q{x,t\y,t') 



-"j exp 



(17) 



where 7-"^ signifies path ordering. Non-eikonal corrections in Eq. ( |16[ ) reflect addi- 
tional Brownian motion along the trajectory. 

Let us now consider the emission of a gluon and for the moment we are only 
interested in the energy spectrum, postponing the discussion of the transverse mo- 
mentum dependence to Sec. |4.4[ The square of the amplitude is represented in Fig. [2] 
(left panel), where the upper lines corresponds to the amplitude and the lower to 
the complex conjugate. The splitting occurs at the time ti and <2 in the ampli- 
tude and the complex conjugate, respectively. After integrating over the transverse 
momenta of the outgoing gluons the regions [0,ti] and [t2,L] cancel out reflecting 
the fact that the production process is fully determined by the three-point function 



between ti and t2^ involving the product of 3 propagators Q, defined in Eq. (16), in 



the transverse plane, see Fig. [2] (right panel). The center-of-mass of the system of 
the three particles, labeled 1 and 2 in the amplitude and in the complex conju- 
gate, is conserved and after integrating over transverse momenta it reduces to zero, 
i.e., zri + (1 — z)r2 — ro = where z and 1 — z are the -I— momentum fraction 
of the parent parton carried by the offspring partons 1 and 2, respectively. As a 
result, the cross-section will only depend on a single transverse separation u, where 
ri — ro = It, ^2 — ro = zu and t"i — r2 = (1 — z)u. The medium-induced gluon 
radiation spectrum reac 



dl 



ind 



a. 



zPR^giz) 



dz 



[z(l — z)p'^ 



Re 



dU 



dti 



X 9„, • 9, 



G{u2,t2;ui,ti\z) - Goiu2,t2;ui,ti\z) 



UI—U2—0 



(18) 



which is fully determined by the three-point function G{u2,t2',ui,ti\z) and where 
Pfj^g(z) is the Altarelli-Parisi splitting function for a parton in a representation R 
converting into a gluon. It is found after taking the appropriate medium averages, 
see Eq. to be given by 

r*2 



G{u2,t2;ui,ti\z) 



Vu exp 



dt 



. z(l — z)p'^ 



1 



M^- -r3(M,t) 



(19) 



April 15, 2013 2:11 WSPC/INSTRUCTION FILE MMT"review'final2 



Jet physics in heavy-ion collisions 15 








z 








\ 1 






'•() 1 



Fig. 2. Left; Square of the gluon emission amplitude in medium. Right: Structure of the squared 
amplitude after integrating out all transverse momenta. 

This three-point function can also be rewritten in terms of Schodinger equation, 
which allows for an expansion in terms of the medium density,'^ 



d 



G{x,t-y,t')^i5{t-t')5{x~y), (20) 



dt 2z{l-z)p+ ' 2' 
where the interaction potential of the 3-body system, see Fig. [2] (right panel), reads 

r3(r, t) = ic^ra (r, t) + (cr - ^Ca^ Fa {zr, t) + ^C^Fs ((1 - z)r, t) , (21) 



propagator collapses to the free propagator in vacuum 



where F2 is given by Eq. (|10|). In the absence of medium effects, F3 — > 0, this 

(22) 



Go{u2,t2;ui,ti\z) 



z(l — z)p^ 



27r(t2 -ii) 



exp 



iz{l — z)p' 



(m2 - 



which has to be subtracted in Eq. ( 18 1 to obtain the purely medium-induced spec- 
trum 

It is worth noting that in the limit of soft gluon emission, i.e. 2: <C 1, the 
three-point function collapses to a two-point function which is only sensitive to 
the broadening of the emitted gluon, since rg — > r2 in Fig. [2j In this case, the 
three-body potential simply picks up the correct color factor for gluon emission, 
FaCr) ~ C^r2(r). 

As discussed in Sec. [3j a more detailed analysis of the interaction of the hard 
parton with a thermal medium can been performed in the framework of the HTL 
approximation. The main quantitative improvement is the correct treatment of 
the medium potential entering in F2 and dressing of the partonic propagators by 



thermal masses. The latter effect can easily be included in Eq (201 by replacing the 



>A more general regularization prescription was introduced in Ref. 41 
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derivatives as follows 

z(l — z)p'^ 



z{\ — z)p'^ 



(1 — z)p+ zp+ p+ 



(23) 



where m.^ is the thermal mass of parton i participating in the emission.^ 



mmzi 



Integrating Eq. (181 is a difhcult task in general and can only be performed 



numericallyjiS^ but one can recover analytical results in some limiting cases. One 
strategy is to expand the three-point function in powers of the medium density, 
see Eq. (20 1. This procedure is commonly called the 'opacity expansion EUEll and is 



applicable in the case of a dilute medium or for short emission times. Whenever the 
number of scattering centers becomes large, one is obliged to resum the contribution 
from all orders in opacity. 

Another strategy is to pick up the leading-logarithmic behavior, as discussed in 
the context of Eqs. (11) and (|9]). Hence, in the 'harmonic oscillator' approximation 
the integrals over the transverse coordinates become Gaussian and can be performed 
analytically. Assuming a medium of constant density, q(t) — qoQ{L — t), one finds 
the well-known BDMPS spectrum 



dl 



ind 



as 



dz 



zPe 



,(z)ln 



(1 + ^1 



2z(l - z)p+ 



(24) 



where Pii^g{z) stands for the relevant Altarelli-Parisi splitting function'"^" for gluon 
emission off a quark or gluon, R — q,g, respectively, and 

q,S = [(1 - z)Ca - z'Cr] qo . (25) 

Solutions for a wide class of expanding media are also analytically available)^li25l 
see also Refs. |109yil0| for a general scaling law. 

For large media or, equivalently, for gluon emissions meeting the condition 
zp+,(l — z)p+ <C we can simplify the spectrum in Eq. (24) further and ob- 
tain an emission rate per unit length dL, which reads 



as z/C(z,p+) — z 



dl 



ind 



dz dL 



2z{l - z)p^ 



-zP, 



iz), 



(26) 



which in this approximation is constant. This reflects the fact that gluons can be 
emitted continuously along the travelled path. Next-to-leading logarithmic correc- 
tions have also been found analytically.'^ This rate is quite useful for phenomeno- 
logical applications, see Sec. |4.3[ 

In the high-energy limit one can neglect the energy of the emitted gluon com- 
pared to that of the initial hard parton. This amounts to putting z <C 1, and one 
typically denotes uj = zp'^ . The spectrum in Eq. ( 24 ) displays two limiting behav- 
iors. When oj <C Wr, 



duj 



%0 ' 



(27) 
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which closely resembles our estimate in Eq. (15 1. It reflects the fact that coherence 
effects build up when the branching time increases causing the suppression of the 
spectrum. When a; 3> Wc, the radiation time tbr has reached its maximal value L 
and the spectrum is even more suppressed, namely 

dr^"^ 2a, C„ -^2 



^ -^P^ — , (28) 
duj 127r V / ' ^ ' 

Going beyond the leading-log approximation and taking into account the correction 

from the short distance behavior of r2 would reduce the power of the suppression 

to ~ (wc/w)^ ESESiiil The leading-log approximation is valid as long as cj > wbh, 

where the soft divergence has to be regularized by hand, see the discussion above. 



4.3. Energy loss 

The mean energy of a hard parton with p"*" ^ ujc 



is obtained by integrating the 



BDMPS spectrum, Eq. (24 1, yielding 



dl'™'^ 

did CO — ; — cx asC„uJc 
duj 



(29) 



up to logs of the energyEHiiiUlii^ It shows the scaling as opposed to the collisional 
energy loss that scales as L, cf. Sec.|3] Thus, for large enough media radiative energy 
loss dominates over collisional one. 

In order to compute the quenching of spectra of high-p^ particles one has to go 
beyond the mean energy losJI^and take into account multiple gluon emissions. The 



rate of emissions, Eq. (26), takes a particularly simple form for soft gluons radiated 
in the medium. In this limit, z <^ 1, we can neglect the requirement of exact energy- 
momentum conservation and thus it is expected to follow a Poisson distribution. 
This assumption also relies on the requirement that multiple gluon emissions off 
the hard parton are independent. While this approximation seems unjustified at 
present, it is actually valid in a specific kinematic range which will be discussed in 
the next section. Then, the probability of emitting a total energy e in course of an 
arbitrary number of emissions can then be written a^HHlHHllll^ 

dP""^ 



P{e) 



ri=0 



1 



n 

.i=l 



dojj 



dujj 



(30) 



where the total number of emitted gluons is iVg = J duj dP'^'^/doj. The limits on the 
w-integration have already been discussed above. In Mellin space, the summation 
can be performed analytically, and a solution found ioP^ 



P{e) 



Zm 



(31) 



where the contour C runs along the imaginary axis and where the Laplace transform 
takes a compact form 



P{i^) = exp 





f , dl'"''^ 






do. 




J 




soft 



(1 



(32) 



April 15, 2013 2:11 WSPC/INSTRUCTION FILE MMT"review'final2 



18 Y. Mehtar-Tani, J. G. Milhano, K. Tywoniuk 



The energy-loss probability defined in Eq. (30) has many applications. As an exam- 



ple, let us illustrate the bias toward small energy losses mentioned above defining 
the so-called quenching factor,'^^^' 



(33) 



{TAB)daPP- 

which compares the production of hadron h in heavy-ion collisions to proton-proton 
collisions at the same energy. Since the vacuum spectra typically drop like a power 
of the transverse momentum ~ p^", the support of the integral in Eq. (33) is biased 
toward small e. Assuming a constant and large slope parameter n, the quenching 
factor is found to scale as Q^{pt) — exp[—AN{pT/TTn)], where AN{lu) is defined 



in Eq. (36). This indicates that only gluons with uj < px/n, which is typically 



smaller than ujc, contribute to a typical radiative process. 

The final distribution of partons with energy E can be evaluated by convoluting 
the distribution P with an initial distribution. 



DiE) 



deP{e)D{E 



(34) 



where the initial condition can, e.g., be taken from pcrturbative QCD, see Sees. [5] 
and |6] for a further discussion. Furthermore, this distribution obeys the following 
integro-differential equatiorHUliill 



dL 



DiE) 



de 



dl 



ind 



de dL 



D{E 



soft 



de 



dl 



ind 



de dL 



DiE) 



(35) 



soft 



which follows by utilizing the Laplace transform in Eq. ( 32 1 . Without going into 



details here, it is possible to demonstrate that as long as the gluon emissions are 
independent this rate equation holds in general for any z, so that one can substitute 



the spectrum in the soft limit by the full rate Eq. ( 26 ) in Eq. ( 35 ) and also include 



mixing between the partonic species HHlii^lll^ In addition to parton splittings, this 
rate equation also includes absorption of thermal particles. In the general case. 



P(e) is not anymore described by the Poisson ansatz as in Eq. (30). In the following 
subsection we will indeed demonstrate that this assumption can be proven in a 
well-controlled way. 



4.4. Factorization of multiple gluon branchings 



From Eq. ( 27 ) one can estimate the average number of gluons emitted with energies 

, dP'"^ _ L 



larger than a given value w, 

ANiuj) 



duj 



duj' 



th-c 



(36) 



where a = asCn/i:. As long as ANiu)) < 1, i.e. w > ci^i^ci it may be identified 
as the probability to emit one gluon with energy lo' > uj. This is the case for the 
relatively hard emissions that dominate the energy loss, see Eq. (29). In this regime 
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the probability of multiple emissions is small. But for sufficiently soft gluons, such 
that 

w < 6?ijJc , (37) 

AN{u!) becomes larger than unity and multiple emissions are clearly important 

In technical terms, this means that all powers of aL/thr have to be resummed. 
By contrast, the soft and coUinear divergences of vacuum radiation give rise to the 
logarithmic enhancement in terms of the jet energy, as log^ E. We ignore this for 
the time being by formally setting L — )■ oo, and postpone a discussion of this issue 
to Sec. E 

The treatment of multiple emission is a priori complicated due to interferences 
between various higher-order processes. For instance, for radiation in vacuum these 
interferences are essential and give rise to angular ordering, see Sec. [2j This mostly 
affects the emissions of soft gluons which typically have long formation times, 
t[ ^ 1/6'^uj. On the contrary, soft gluons induced by the medium are produced 



very rapidly, see Eq. (13). In fact, for the regime of multiple emissions, defined by 
Eq. (37), we naturally obtain that thi L — allowing us to treat these emissions 
as quasi-instantaneous. Below, we will also show that thr represents the time when 
the products of the branching decohere in color ,1^211 see Eq. ([ZTjlp] In other words, 
after a time ~ ^br the partons formed in the branching propagate independently of 
each other and all interferences are suppressed, see Sec. [5j 

Consider two successive branchings. The smallness of tbr severely limits the 
available phase space for interferences, which scale like AA^™* oc {asL){astb^) due 
to the requirement of overlap close to the branching. The factorizable piece, on the 
contrary, scales with the maximal phase space, ATV^^'^' c>c (agL)'^, due to the quasi- 
instantaneous nature of the radiation. In line with the discussion above, it follows 
that the interferences are suppressed as compared to the factorizable contribution 
when ^ by the factor 

AA'"t ibr 

In order to address the issue of multiple branchings in more detail in the regime 
discussed above we shall proceed by analyzing the two first orders in a parton 
cascade which will allow for a generalization to all orders. We shall also extend the 
discussion of the previous section by discussing quantities that are fully differential 
in transverse momenta. 

Consider a gluon produced at an initial time by a hard process described by 
the cross section dahmd/d^po' where dilk = (27r)~^(i^pgdpJ"/2p(j" is the standard 
phase space measure. The cross section for observing the gluon at some later time 
with transverse momentum k is obtained by convoluting the initial spectrum at 
time to with the probability for propagation from to to t^ which is described by a 



''For a more general discussion of this physics, we refer to Sec. [H] see in particular Eq. 



49 



April 15, 2013 2:11 WSPC/INSTRUCTION FILE MMT"review'final2 



20 Y. Mehtar-Tani, J. G. Milhano, K. Tywoniuk 



similar evolution equation as Eq. (20) with Fa replaced by CAr2- The result reads 



d(7o f <fiPo , , ^ rfChard 



^rik~Po;tL,to) ^j^, (38) 



where V{k—pQ; ^l, io) represents the probabihty that the gluon acquires a transverse 
momentum q = k — from the medium in course of the propagation. This is the 
first building block of our probabilistic picture for jet evolution. In general it is given 

by 



V{q;tL,to)^ Jd^rexp -iq ■ r - ^Ca dt'T2{r,t') 



(39) 



For homogeneous media and in the "harmonic oscillator" approximation, valid as 
long as |<7| ^ q(iL — to), it takes a particularly simple form 

4^ f q^\ 



^(.,At). — exp^-^j, (40) 

where At = <l — to- In the opposite regime of 'hard' interactions with the medium, 
see discussion around Eq. (Ill, this distribution behaves as V{q,At) ^ q This 



classical propagation in the medium is characterized by the global scale — qL, 
corresponding to he maximum transverse momentum that a particle can acquire 
since typically At ^ L. Thus, due to the broadening in medium the propagating 
particle tend to deflect to large angles. See also Sec. [3] for a further discussion of the 
role of q. 

We now turn to the cross section for the splitting of an initial gluon into two 
offspring gluons, labelled a and b. This process is depicted in Fig. We distinguish 



three time intervals: (i) the propagation before the splitting, described by Eq. (38 1 



(ii) the splitting process itself, described by a three-point function, as discussed in 



Sec. 4.2 and (iii) the subsequent propagation of the produced pair. 



Analogously to Sec. 4.2 it can be shown that emission region (ii) has support 

at most of the order of the branching time tbr = \J z{l — z)pq /qcS^^ Since we are 
working in the limit ibr ^ this region is replaced by a quasi-local vertex in Fig. [3] 
describing the quantum branching process. 

After the products of the branching are formed, their further propagation in re- 
gion (iii) involves their correlations in color and momentum space and is in general 
described by an, a priori unknown, four-point function ^^^^(ti,^). The key result 
of the factorization proof,'^^ is to demonstrate that the four-point function factor- 
izes into independent propagations of the two gluons in the limit tbr ^ i up to 
corrections of the order of ibr/^ 



^(4) {t^^t) (X V{ka -p,tL~ t)V{kk -q + p,tL-t)+0 (^-f^ 



(41) 



This proves the decoherence of induced gluons shortly after they are produced. 
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Fig. 3. Graphical illustration of the equation |42| . The thick wavy lines represent the probability 
V for transverse momentum broadening, the black dot is the splitting probability K, and the 
circled cross is the cross section of the hard process producing a gluon of momentum po • 



i.e. at times larger than tbr- The resulting spectrum at leading order in i/ibr reads 
theiPSl 



— = 2g^z{l - z) 



X lC{p- zq,z,p+,t)V{q-pQ,t-ta) 



At I V{ka-p,tL-t)V{kb-q + p,tL-t) 

PO'I'P 



da 



(42) 



Po 



where z = k'^ /pq and we have adapted the notation Jp = J cZ^p/(27r)^. This result 
can be interpreted as a classical branching process, as illustrated in Fig. [3} in the 
following sense. After propagating from to to t, during which it acquires a transverse 
momentum q—pQ, the original gluon splits into offsprings a and b with a probability 
~ aslC{p— zq, z, g+) which depends on the longitudinal momentum of the parent 
parton, the longitudinal momentum fraction z = p'^ /q^ carried by gluon a, and the 
transverse momentum difference p — zqQ After the splitting, the two gluons a and b 
continue to propagate through the medium, from i to i^, thus acquiring additional 
transverse momentum. 

The quasi-instantaneous, fcj_-differential splitting kernel can be computed, sim- 
ilarly as in Sec. |4.2( yielding 



}(^{p,z,Po,t) 



[zil-z)p^ 



Re 



dAt e-'"^ P 



■ [G{u2,t + At;ui,t, z) - G'o(m2,< + At;ui,t, z)]^^^^ , (43) 



cf. Eq. (18). For a homogeneous medium the splitting kernel is independent on the 



emission time t. One can go further an evaluate Eq. (43 ) in the "harmonic oscillator" 
approximation, where it read^^^H 



l^iP,z,Po 



2(1 - z)Pq 



■Pgg(z) sin 



PL 



exp 



PL 

'^ki- 



rn 



'The conservation of longitudinal momentum implies of course pg" = 5+ = fc^ 
p+ = 2q+. 



+ k+ with k2 
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where k^^ — \J z{l — z)p(|goff is the typical transverse momentum generated during 
the branching process. This kernel generalizes the result obtained in Ref. |117| in the 
eikonal limit. 



The spectrum in Eq. (42) established the factorization of induced gluon radi- 



ation in medium and constitutes the building block for jet evolution in the limit 



th-c <C i at large angles, away from the jet core. Note that integrating Eq. (44) over 



the transverse momentum we recover the branching rate defined in Eq. (26). Hence, 
the independent branching approximation assumed in the rate equation introduced 
irl25lliiil is justified so long as the typical in-medium branching time is much smaller 
than the medium length and can be generalized by including the transverse mo- 
mentum dependence of the in-medium shower.l^^ 



5. Color coherence and color flow in the medium 

Having discussed the main elements of the physics of jet evolution in vac- 
uum, on one hand, and the dynamics of medium-induced radiation, on the 
other, we can now combine these insights to study how the intrinsic jet- 
structure is modified by the QGP. For the sake of simplicity, these aspects 
will mostly be discussed in the context of the smallest building block of the 
jet, the antenna system, where one most easily can illustrate the effects of 
color decoherence. 

So far we have discussed aspects of the branching of a highly energetic particle 
separately in vacuum and in medium (where in the former case the particle must 
be strongly off-shell to be able to emit). The complications arising from resum- 
ming multiple gluon emissions were treated in two completely different manners in 
the two cases. In the former case, the probabilistic interpretation of the shower at 
leading-logarithmic approximation is saved by the condition of angular ordering of 
subsequent radiation, see Sec. [2j In the latter, on the other hand, at least in the 
regime of soft gluon radiation, i.e. uj < u^, one can rigorously prove the complete 
factorization of subsequent radiation in the medium, see Sec. |4.4| Thus, while a 
vacuum jet is a completely color coherent object the shower created via medium- 
induced branchings is completely decoherent. The natural evolution variable of the 
two types of shower is also different: in the vacuum the evolution is ordered in vir- 
tualityp^ see Sec. [2j while in the medium the emission rate of independent gluons 



scales with the (physical) length of the medium, Eq. (26). 

Dealing with jet systems with an intrinsic scale much larger than any scale re- 
lated to the medium, e.g. E = px 300 GeV/c at LHC, we should therefore proceed 
to discuss these seemingly disparate physical processes in a common framework. A 
complete theory of jets in heavy ion collisions would be dictated by the interplay 
of the relevant hard scales in the problem. While the dynamics of jets in vacuum 



Equivalently, in transverse momentum or angle. 
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is fully determined by the transverse mass of the jet Q = EQ, where Q is the jet 
opening angle, and the non-perturbative scale Qo, in the presence of a medium of 
length L two additional global scales arise: the typical transverse momentum ac- 
quired by a particle traversing the medium Qs = y/qL, which inverse relates to the 
resolution power of the medium to color charges in the transverse plan, and the 
typical inverse size of the jet in the medium — (8i)^^. As we shall discuss 
in what follows, the competition between these two scales determines the degree of 
alteration of color coherence of the parton shower. 

Naturally, in addition to inducing radiation, the presence of a deconfined medium 
— i.e., a medium interacting with the probes via color exchanges — is also expected 
to alter the coherent structure of the jets. Thus, the notion of a resolution scale of the 
medium needs to be established. This was first studied in the context of radiation 
off a color correlated pair of partons, the so-called antenna configuration,'^!!!!!!] 
traversing the QGP|3 

For two color-correlated probes propagating through the medium with momenta 
Pi and p2 , the probability to lose their color-correlation after a certain (light-cone) 
time t prior to further gluon radiation, is described bjliillEillliHH] 



Amcd(0 = 1 - exp 



dt'T2{5n^2t',t') 



(45) 



where |(5ni2| ^ 
is given by Eq. 

so that the second term in Eq. ( 16 ) can be neglected. In the 'harmonic oscillator' 



9 12 is the angle between momenta of the constituents and r2(r,t) 
[TOl We have assumed that the particles are both energetic enough 



approximation this expression takes a particularly simple form, involving only the 
characteristic hard scales of the problem, namely 



A„icd(t) = 1 - cxp 



(46) 



where rj_{t) — 6i2t and Q'^{t) — qt. The medium characteristics are all encoded in 
the parameter q, discussed at length in Sec. |3] The function in Eq. (46) behaves as 
an order parameter describing two characteristic regimes and is therefore referred 
to as a 'decoherence parameter'. In the two extreme limits, it reads 



A„,ed(i,0) 



, 7-j_ <C ^ (coherence) , 
1 , 7'_L ^ Qj^ (decoherence) , 



(47) 



which possess a simple geometrical interpretation, see Fig. |4] In the former case, 
sometimes called the 'dipole' regime, the pair remains correlated after propagat- 
ing for a time t. In this case one can expand the exponent in Eq. (|46|) to obtain 



A 



mcd 



which is the characteristic behavior of color transparency^ 



the latter case, called the 'decoherence' regime, the pair is resolved by the medium 



''Note that the color charge of the pair is not necessarily balanced, so that the antenna can possess 
a total charge. 
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Fig. 4. The characteristic regimes of radiation in media: the 'dipolc' regime, r± <g; Qs ^ (left) 
and the 'decoherence' regime, r± 2> Q7^ (right)- Figure taken from Ref. |123| 



interactions. In particular, this implies that all interferences between the particles 
are suppressed and the two particles behave as independent color charges. 

The role of color decoherence is most easily discerned in the soft sector where 
medium-induced gluon radiation is suppressed. Let us, for the purpose of illustra- 
tion, presently focus on this regime. In this limit, all medium effects enter via a 



multiplicative factor Amod , defined in Eq. ( 46 ) , and the coherent spectrum off one 



of the antenna constituents, i.e. the spectrum in the presence of the charged carried 
by the other constituentj^ after integrating over the azimuthal angle, can be written 
as 



agCji did d6 



ei 



^mcd 



ei 



2) 



(48) 



where uj and 9 are the energy and angle of the emitted gluon, respectively. These 
features are illustrated in Fig. [5] In the absence of a medium, Amed 0, we recover 
the pure vacuum spectrum, which is explicitly constrained to angles smaller than 
the opening angle of the pair. This reflects the property of angular ordering.^^ 
For a finite medium density, Amcd > 0, the medium-driven component occurs at 
large angles and is suppressed inside the antenna.'^^ In the limit of a completely 
opaque system, Amed — ^ 1, the interferences are completely washed out and the 
soft emissions in the presence of a medium reduces to independent radiation off the 
quark and antiquark, as if they were radiating in the vacuum. From this simple 
picture we see that one of the chief implications of the onset of decoherence is that 
it marks the breakdown of the characteristic angular-ordering features of the jet as 
it propagates through the medium. 

Correspondingly, the same condition of decoherence determines how the medium 
'resolves' the effective charge that subsequently will loose energy via induced radi- 
ation, described in detail in Sec. |4|iiZlliH^ For higher gluon energies, following the 
above discussion, the interference diagrams also involve contributions from purely 
medium-induced gluons. The emerging picture can be summarized in the following 
simple principle. 



°For a definition see Ref. |T24|and|117][l23| 
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Fig. 5. The soft gluon emission spectrum off the quark constituent of a singlet antenna with 
opening angle 0\2 = 0.2, according to Eq. ( |48[ l, in the presence of a medium with Aj^od = 0.5 
(solid line). Here a = ctsCp/ir. On average vacuum radiation is confined within 6 < 9i2, while the 
medium-induced radiation is radiated at 6 > 612- The limit of opaque medium, given by Amcd = Ij 
is marked by the dashed line. Figure taken from Ref. |119| 



Amcd = — > medium-induced radiation off the total charge, 

Amod = 1 ^ medium-induced radiation off the constituent charges. 

For instance, in case of an antenna in a color singlet configuration, e.g. originat- 
ing from the splitting of a virtual photon 7* — qq, no radiation is induced by 
the medium in the 'coherent' regime. Consequently, we observe that the decoher- 
ence of the color correlation between partons propagating through the medium also 
determines how they loose energy via radiative mechanismsj^ 

The dynamics encoded in the decoherence parameter can also be translated in 
terms of a characteristic decoherence time, 

where in the last step we multiplied both numerator and denominator by the the 
energy carried by the antenna constituents. The timescales t[ ^ 1/E0^2 ^-nd fbr, 



see Eq. (13), are the familiar branching times in vacuum and medium, respectively. 
Here, they both relate to the timescales over which the antenna itself is formed. This 
point brings us to two crucial insights. For highly energetic, formally — > 00, and 
strongly coUimated probes, 612 <^ 1, the timescales for formation and decoherence 
are strongly separated, i.e., t{ <ti t^. Such configurations can only be realized by 
vacuum emissions. Secondly, for radiation which is characterized by tf — tbr, we 
observe that if = t^. Since this type of radiation typically is distributed at large 
angles away from the jet direction, see Sec. [4| we will not discuss it further in this 
Section. 



PNote that the same coherence effects are at play for purely elastic interactions in the medium.li^SI 
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Since the description of jet fragmentation is governed by the relevant scale, or 
ordering variable, it will be helpful to rephrase the physics of decoherence in the 
same language. Let us for the time being consider global scales, i.e., involving the 
total length of the mediumj^] Thus, besides the running scale of the intrinsic jet 
evolution, in the presence of a medium a novel hard scale given by 

Qmod = max (rJ^^jQs) , (50) 

which determines the range (in transverse momentum) of the influence of the 
medium, will come into play. In particular, radiation with k± > Qmcd is indeed 
coherent. This condition is met in two distinct cases. 

The first case concerns large-angle radiation. If we consider the ordering features 
of radiation at angles 6 > 9i2, the discussion above describes how angular ordering 



is spoiled for soft gluon radiation by the medium decoherence effects, see Eq. (48 1. 
For finite gluon energies, the maximal angle of radiation off the constituents of the 
pair becomes Qmed / 11111112^ This angle replaces essentially the opening angle of the 
pair for subsequent emissions. In the 'decoherent' regime, where Qmcd = Qsi this 
angle agrees with the maximal angle of medium-induced radiation which provides 
a subtle consistency check: the system cannot transport its color charge to larger 
angles. 

In a complementary fashion, coherence effects are also important for very 
collinear radiation. As pointed out above, these configurations are characterized 
by tf <C id- Due to the ordering features of the vacuum shower, this opens up the 
possibility for multiple emissions within the window delimited by t^ which remain 
color coherent. This applies generally to hard and collinear radiation inside a jet, 
specifically all modes of the jet which meet the criterion for coherence .1^221 This possi- 
bility is particularly appealing in the context of high-p^ jets. The presence of a hard 
scale due to the medium therefore converts to the following angular condition. All 
jet sub-structure contained within the critical angle which, in the approximations 
used throughout, reads 

e,= {qL^)-"\ (51) 

remains unresolved by the medium. Analogously to the antenna case, resolved 
jet substructures, i.e. separated at angles larger than 9c^ evolve independent of 
each other and are not anymore forming a color coherent structure. Of course, 
this effective number of resolved subjets also depends on the specific fragmentation 
property of a given jet. It follows that the jet evolution on angles smaller than Oc is 
completely vacuum-like and, consequently, one should expect that the distribution 
of these fragments scales like the vacuum distribution. Secondly, the unresolved 



system looses energy coherently with a rate as in Eq. ( 26 1 proportional to the charge 
contained within Oc- 



iWe will denote = r^{L) = 012^ and Qs = QsiL) = 
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For realistic situations it was notecP^ that a core containing most of the jet 
energy remains unresolved by the medium. This makes it a particularly appealing 
probe of the transport parameter q of the medium. Secondary substructures carry 
typically smaller energy fractions and are expected to be influenced quite violently 
by medium effects. In addition to the effects on final-state (time-like) showers, color 
decoherence is also paramount when studying processes which involve several color 
charges, e.g., observables which become important at NLO, see Sec. [6] for a further 
discussion. One particular example are multi-jet events, which in hadronic collisions 
typically involves the interplay of space-like and time-like emissions. ^■^^ An in- 
teresting prospect would be to observe the modification of their coherent structure 
in the presence of a medium .1^221 

As discussed above, decoherence in course of subsequent branching leads to the 
disturbance of the color flow within the jet. This influences the later stage of the 
showering and affects the particle distributions created in the course of hadroniza- 
tion.li^SI We emphasize that the established models of this non-perturbative process 
contain the basic elements of color coherence. That interactions of an outgoing probe 
with soft gluons from the background 'reorganize' the color flow of the event has 
also been realized in the context of diffractive deep inelastic scattering |i3i|li32j j^g_ 
cently, it was noted that constituents of a jet that have interacted with the medium 
could end up in a configuration where they are color-connected with the medium 
excitations, rather than to other jet fragments. ^^^Hl These modifications have been 
calculated up to second order in medium density and also checked numerically with 
existing event generators.'^^ This effect causes the formation of strings or clusters 
with typically larger invariant mass than for unperturbed color flow as in vacuum, 
leading to larger multiplicity of soft particles and consequently to overall softer 
spectra and a characteristic high-p^ depletion. 

6. Phenomenology of jet quenching 

We have reviewed how decoherence of the jet governs its fragmentation pat- 
tern in the medium and determines the rate of coherent energy loss in some 
limiting cases. Based on these insights, we try to assess the present status of 
jet quenching phenomenology and the challenges ahead. The wealth of high- 
quality experimental data can help in pinning down the relevant effects. 

6.1. Modeling the shower 

The focus of this review is the physics of high-pr probes, and jets in particular, 
which traverse the complex state created in the wake of heavy-ion collisions. So far 
we have reviewed some of the chief elements related to parton branching, such as 
accumulation of transverse momentum and color decoherence, but the theoretical 
framework unifying these various aspects and allowing for a well-controlled com- 
parison with data is yet to be established. A key issue currently under investigation 
is the space-time picture of the branching process and, in particular, how collinear 
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and medium-induced radiation intertwines in the shower. For phenomenological ap- 
phcations one currently employs working models which are based on reasonable 
assumptions within a certain window of kinematical validity. Considering the lim- 
ited space available, here we will only comment on a few key issues. 

The single-gluon emission spectrum in a static medium defines a clear theo- 
retical setup — often referred to as the 'QCD brick problem' for jet quenching.Esl 
The main input to this calculation is the medium interaction potential containing 
the thermal effects and the mean free path, discussed in Sec. [3] In the high-energy 
limit this allows to compute, strictly speaking, the energy spectrum of medium- 
induced gluons, albeit numerically.'^^ Instead, most model implementations employ 
well-controlled approximations, such as the 'harmonic oscillator' approximation dis- 
cussed extensively above, which allow for analytical solutions. In this case, care has 
to be taken to correctly incorporate the limits of applicability and estimate errors. 



For instance, for a fixed value of the quenching factor, defined in Eq. (33), various 
implementations report a wide discrepancy of the resulting medium parameters, 
mainly q.^^ ^^'^ This reflects the theoretical uncertainties inherent in the employed 
approximations]^ These uncertainties should be systematically improved upon. In 
particular, finite-energy corrections can also play a significant role for quantitative 
estimates |92|[i35|[T36| j^ecently, a numerical Monte Carlo procedure has been formu- 
lated to correctly reproduce the LPM effect in QCD.^i^ 

Due to the large phase space for radiation for high-energy probes, the second 
great challenge is to implement multiple emissions in the medium. It is useful to 
group the current working models into two broad categories set apart by the choice 
of evolution variable for the shower. 



On one hand, the constant medium-induced rate, Eq. (26) and the discussion 



below, seems to indicate a strict time ordering of medium-induced radiation. The 



branching in medium is then described using the rate equation Eq. ( 35 ) , generalized 
to account for energy-momentum conservation and including mixing of partonic 
species. Implementations of this scenario usually assume that the input partonic 
distributions can be taken as vacuum HHliissj 

On the other hand, assuming that the initial virtuality of the jet is much larger 
than any of the scales from the medium, one can treat the medium-induced split- 
ting processes as a small perturbation on top of the vacuum evolution ordered in 
virtuality which resums coUinear emissions, see Sec.[2j This approach is also appeal- 
ing from the point of view of Monte-Carlo implement at ion CMMIj ^^^^ argue 
heuristically that this holds for splittings where the drop in virtuality exceeds the 
amount of transverse momentum accumulated from the medium during the for- 
mation time,'^ k'j_ > Vq^j where k± is the transverse momentum of the gluon 
created in the splitting. This can be accommodated in various ways. Here we bring 
forward one particular prescription where a medium-modified splitting kernel is 



■^Similar studies, including an expanding medium, have also been documented in Refs. |93||94| 
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definedHmiEllElll 

Pl°\z) = + ^P^^\z; . . .) , (52) 

which generahze the standard Altarelh-Parisi^° ones. The medium-modified part, 



second term in Eq. (52 1, is again related to the rate of medimn-induced giuons. 



see e.g. Eq (24) and Refs. 43 44 144 for other deriviations. Note that it typically 
involves the global scales involving the full medium length L, setting it apart from 
the rate equations discussed earlier on, which are local in time. The full splitting 
function is used as input to the standard DGLAP evolution equations of the hadron 
fragmentation function, see Sec. [2] In the limit of soft giuons, z <C 1, this approach 
also leads to nearly independent emissions .1^^ 

The elements discussed so far can be improved systematically, at least in some 
limiting cases, within the framework of perturbation theory. Additionally, the in- 
fluence of several non-perturbative effects are necessary for a realistic comparison 
with data and the complexity of the problem at hand grows rapidly. This purely 
phenomenological input, from the point of view of jet quenching, includes the mod- 
eling of the initial nuclear geometry and space-time picture of the plasma evolution, 
typically taken from hydrodynamical models. Considering jet quenching as a tomo- 
graphic probe of the QGP, these effects can play an important role when comparing 
to experimental observables.'^^ The fate of the produced color charges in course of 
medium-induced radiatiorP^ and their back-reaction on the plasma evolution is still 
under discussion. Finally, as discussed in Sec. [5j medium effects on hadronization 
can also influence the measured hadronic spectra.'^^ An additional complication 
arises for jet observables which demand subtle background subtraction procedures 
for an apples-to- apples comparison to experimental data.'^^^HHH 

Presently, state-of-the-art calculations of jet quenching effects usually involve 
numerical event generationilUlH^I ^'^^ll^'^^ lll^ll^ Although many codes are used by 
several groups for comparison with experimental data, only a few are properly 
documented in detail. 

6.2. Lessons from the data 

With the advent of the heavy-ion program at the LHC the physics of heavy-ion 
collisions has entered a new, high-energy regime]^ The unprecedented kinematical 
reach together with a large lever in collision energy compared to RHIC provides a 
wide range to test the physics of energy loss under widely different conditions. Ample 
production of high-p^ jets and the excellent detector capabilities allow for the flrst 
time precise experimental data on their modifications in various aspects. Without 
going into the details of jet reconstruction in heavy-ion collisions the jet 
sample discussed below usually comprise high-pT jets reconstructed using the anti- 
kt algorithniUllESIl with a relatively narrow cone size R ^ 0.2-0.5 to minimize the 



^Since the RHIC heavy-ion results have already been reviewed elsewhere, see e.g. Ref. |15[ we will 
presently focus on the now results from the ALICE, ATLAS and CMS Collaborations at the LHC. 
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effects of background fluctuations. For an overeview of the most recent developments 



together with references to unpublished results, see also Ref. 156 

Just as for the inclusive hadron production, a key measurement is the inclusive 
spectrum of jets in heavy-ion collision. Due to the highly nontrivial procedure to 
extract the jets from a violently fluctuating background these unfolded spectra 
are not yet published but many of the inherent uncertainties cancel out in the 
construction of the Raa modification factor Eq. ([T]) (or the equivalent Rcp where 
central collisions are compared to peripheral ones). The measurements indicates 
robustly a suppression by a factor of two for jets with pr = 50-200 GeV/c?^ 
which is roughly consistent with the measurement of the suppression of high-p^ 
charged hadrons |i591|i60| The suppression at the highest available pr is remarkably 
independent of the cone size except at px = 50-100 GeV/c where a decreasing 
trend of the suppression is seen as one opens the jet cone.l^^ These data, together 
with the path length dependence as, e.g., measured by the second azimuthal Fourier 
coefficient V2 of high-pT jets, see Ref. |156[ could help in constraining the nature of 
energy loss for these energetic probes. 

Another important observable is the di-jet energy imbalance.'SSHlII to the 
surface bias of triggered probes one typically expects the sub-leading jet to traverse 
a significantly longer path in the medium than the leading one. Remarkably, while 
the di-jet data manifest a significant enhancement of energy imbalance, the system 
remains almost ideally back-to-back. In fact, in central collisions the sheer number 
of candidate di-jets that pass the required experimental cuts is strongly reduced. 
This would imply that energy is flowing out of the jet cones mostly as soft quanta, 
Pt < 2 GeV/c, and this expectation is indeed confirmed by looking at very large 
angles away from to the di-jet axisESEU This dynamics put strong constraints on 
the mechanism of energy loss — in particular, the mechanism behind the transfer 
of a significant energy into soft modes at large angles. These general features have 
also been confirmed using photon-jet correlations.!^ 

The recent data also allows to study the jet substructure both in terms of lon- 
gitudinal momentum fraction z of jet constituents and in terms of their the angle 
r away from the jet axis per observed jet,!^ thus providing a 'tomography' of jet 
modifications. Due to the QCD dynamics one naturally expects the hardest frag- 
ments to be distributed narrowly around the jet axis while the soft components 
typically tend to extend up to large angles. Thus, these combined measurements 
confirm that the 'core' of the jet, i.e. hard components close to the axis, are not 
strongly affected by the medium. This is generically anticipated due to the physics 
of coherence, see Sec. [5] The soft components at the fringes of the jet, on the other 
hand, are enhanced compared to the vacuum expectation.l^S 

The features of the experimental data follow generic expectations of radiative 
processes in the medium — on one hand, the physics of coherence related to highly 
coUimated and energetic substructures in jets and, on the other, the large-angle 
medium-induced spectrum which carry the imprint of the local medium scale. 
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